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Abstract —Random access is a multiple access communication 
protocol where the users simultaneously communicate with a 
hase station (BS) in an uncoordinated fashion. In this work, we 
consider the problem of multiuser detection in a random access 
bandwidth request context. We propose an enhanced random 
access scheme where the flxed/low-mobility M2M devices pre¬ 
equalize their random access codes using the estimated frequency 
response of the slowly-varying wireless channel. Consequently, we 
have developed two different multiuser detection algorithms. The 
first algorithm works in a greedy fashion where it performs cross- 
correlation of the received signal with a set of decoder sequences 
and detects active users based on the correlation output. We 
derive the condition under which the algorithm can detect a 
given number of active users with high probability. Subsequently, 
we demonstrate an efficient decoder design procedure which 
enhances the user detection performance. A basis mismatched 
sparse recovery technique has been applied in the second 
algorithm which exploit an inherent structure of the random 
access protocol. The performance of the proposed schemes is 
demonstrated in a WiMAX network environment. 

Index Terms —WiMAX, M2M, Sparse representation, Random 
Access. 

I. Introduction 

The orthogonal frequency-division multiple access 
(OFDMA) scheme has been adopted by the IEEE 802.16 
WiMAX standards. The future wireless communication 
network has to support a large number of fixed/low-mobility 
machine-to-machine (M2M) devices that will transmit bursty, 
small data packets (e.g. meter readings, sensor reports etc.) 
under a valid security association with the network HI. This 
corresponds to a heavily uplink-based traffic model following 
a Poisson distribution. Random access based bandwidth 
request (BR) is preferable for such traffic that exploits the 
benefits of statistical multiplexing to support a large number 
of devices with a fixed overhead. In the IEEE WiMAX 
standards, the BR procedure starts with the allocation of a 
predefined set of subcarriers by the BS in a specific time 
slots. The subscriber stations (SSs) which wish to request for 
bandwidth to the BS can take the opportunity by modulating 
a randomly selected code onto the allocated subcarriers. At 
the receiver end, the BS is required to detect the multiple 
active subscribers. Once the BR subscribers are resolved, the 
BS will allocate bandwidth for the corresponding subscribers. 
However, the process will be efficient when the base station 
can separate the BR subscribers successfully. 

The idea of random access channel is related to the concept 
of multiple access channel (MAC) in network information 


theory m. The concept of MAC has been applied success¬ 
fully in CDMA systems 0. However, the application of 
classic MAC channel analysis in the present scenario is not 
straightforward (see ||4l and references therein). The problem 
of multiuser detection (MUD) in a random access on-off 
channel has been studied in a-ia. The work in Q solves 
the MUD problem using a sparse signal recovery framework 
and derives necessary conditions under which the orthogonal 
matching pursuit and Lasso algorithms can detect active users 
successfully. However, the proposed method assumes perfect 
synchronization among all users, which is difficult to guarantee 
in practice. A modified sparse representation framework has 
been proposed in 0 where the above assumption has been re¬ 
laxed. A reduced dimension MUD method has been proposed 
in 0. It has been shown that the computational complexity 
of the reduced dimension MUD is lower than conventional 
correlation based MUD. However, all the methods assume 
that the transmitted signal is subject to flat fading channel. In 
a metropolitan wireless network the communication channel 
generally subject to frequency selective fading and hence the 
above assumption may not hold in practice. Eurthermore, the 
sparse representation of the received signal is not straightfor¬ 
ward for selective fading channel. 

Multiuser detection in IEEE 802.16 based BR has been 
studied in IfTI- lfTOl . A simple solution to enhance MUD is 
to increase the number of BR channels to accommodate more 
users/devices. However, under the existing schemes, only a 
handful of codes can be detected reliably per channel per 
frame in presence of multiple access interference (MAI) from 
different codes, and random noise and frequency-selective 
fading in the multipath wireless channel Q. Therefore, a 
lot of random access channels are required, which would 
substantially reduce the payload capacity of the overall system. 
Hence, random access has been considered as one of the key 
bottlenecks by the IEEE 802.16p working group on M2M 
communications. To overcome this limitation, the recent IEEE 
802.16p amendment has proposed several solutions, mostly 
based on access control over the MAC layer 0. On the 
other hand, a number of works have already appeared in the 
literature concerning this problem, e.g. a and cni. However, 
these schemes require significant modification to the existing 
standards and are not suitable for a network supporting both 
M2M and non-M2M traffic. 

In this work we propose an enhanced random access MUD 
in M2M communication environment. The MUD problem has 
been resolved by two different algorithms. The contributions 
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of the work can be summarized as follows; 

1) We propose a pre-equalized random access scheme. We 
assume that the BR subscribers only have an approxi¬ 
mate knowledge about their channel frequency response. 
The BR subscribers pre-equalize the random access 
codes by using the estimated channel frequency response 
and transmit to the BS. The imperfect knowledge of 
channel frequency response results an additional noise 
term in the received signal at the BS. We analyze the 
statistical property of the noise term and develop a 
suitable data model of the received signal for MUD. 

2) We apply a very simple correlation based approach 
called CMUD for resolving the MUD problem where 
the received signal is correlated with a set of decoder 
sequences. We develop the necessary condition under 
which the CMUD can resolve a given number of active 
users from the received signal. 

3) The necessary condition will show that one can enhance 
the performance of CMUD by controlling some prop¬ 
erties of the decoder sequences. We then propose an 
algorithm to design efficient decoder sequences. 

4) The computational complexity of CMUD is very low, 
but it cannot detect a large number of active users. To 
enhance the user detection performance, we resolve the 
MUD problem by using a basis mismatched sparse sig¬ 
nal recovery algorithm. We exploit some inherent prop¬ 
erties of MUD and formulate an optimization problem 
which can enhance the performance of the underlying 
basis mismatched sparse recovery algorithm. However, 
the optimization problem is non-convex in general. We 
apply a Lagrangian Dual Relaxation method to solve the 
optimization problem. 

Notations: Superscript T denotes matrix transpose. E(a;) 
denotes expected value of x. For a complex number x, its 
real and imaginary parts will be denoted by \x\r and [a:]^ 
respectively. A component of the matrix C at its Uth row and 
j-th column will be indicated by Cij and the j* column of 
C will be represented by Cj. The cardinality of a set T will 
be denoted by ^T. 1^ denotes a vector of length L whose 
all components are one. The fp norm of a vector is defined 
as ||a:||p = \xt\PY^^. diag(a;) refers to a diagonal matrix 
with vector x on its diagonal. 

11. Data Model and Problem Statement 
A. General System Model 

Consider a single-cell WiMAX network with time division 
duplex (TDD) OFDMA physical layer. A BR channel is 
comprised of L randomly chosen subcarriers over one OFDM 
uplink symbol. Suppose the indices of BR subcarriers are 
{jm '■ tn = 1,2,...,L}. When a subscriber station (SS) 
intends to send a BR to the base station then it selects an 
available uplink BR slot and sends a BR packet to the BS. The 
packet consists of a BR code, which is a L-bit pseudo random 
binary sequence (PRBS) chosen with equal probability from a 
bank of K codes HD, where L < K. For rest of the sequel, we 
denote the BR code-matrix as C G where every column 

of C represents an independent code. Each Cij is modulated 


by binary phase shift keying (BPSK) i.e. Cij € {—1,-|-1}. 
The code-matrix is known to BS and every subscriber. 


B. Pre-Equalization 

In this proposed random access model, the SS does not 
send the random access code directly to the BS, instead 
it transmits a pre-equalized version of the code. According 
to the IEEE 802.16 standards, the first OEDM symbol of 
each WiMAX frame is a preamble transmitted by the BS, 
where the subcarriers are BPSK modulated with a boosted 
pilot sequence HD- Typically, the SSs use this information to 
estimate the channel frequency response (CER) for the OEDM 
demodulation process. A SS can pre-equalize its BR code 
using this estimated CER exploiting the channel reciprocity of 
the TDD system EI.a number of pre-equalization techniques 
are available. Eor more details, please refer to HD. 

Consider a time instant when M number of SSs are simul¬ 
taneously contending on the same BR channel. Let, the m* 
SS selects the fcm* column of the code matrix C G 
where m = 1,2,...,M. Considering zero-forcing (ZE) pre¬ 
equalization HD, the transmitted code over the subcarrier 
from the m* SS be 


= forVZG{l,2,...,L}, (1) 

where hi^m is the pilot-aided CER of jf* subcarrier estimated 
by the m* SS. To be more precise, hi^m can be expressed as 

hl^rn — hl^rn “t” (2) 

where /i; is the actual frequency response of the j)* subcar¬ 
rier and ei^m ^ CAf{0, crl is a zero-mean complex Gaussian 
noise with variance 

In the BS, after down-conversion to baseband and OEDM 
demodulation, the received signal from the m* SS over the 
j;* subcarrier be 


yUrn=Xl^rnhl^rn = C>i^kJY^‘^ for V / € {1, 2, L}, (3) 

where is the effective channel experienced by the BS 
due to transmission form m-th SS over the j;-th subcarrier. 
Note that, although the position of the BS and the SS will 
remain almost fixed for stationary/slow-moving M2M devices, 
the channel is continually affected by the movement of the 
external scatterers in the surrounding environment IH . Conse¬ 
quently, will differ from hi^m- To generalize, the effective 
channel can be modelled as na 

^/,m — “f (4) 


where ry/ „j ^ CAf{0, is a zero-mean complex Gaussian 

noise and a is some deterministic complex valued constant. 

Considering the above phenomena, the combined received 
signal at the BS from all M stations over the j/* subcarrier 
will be 


yi 



hl^rn “t” ^l,m 


+ -di, 


(5) 
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where "di ~ CAf{0,a'^) is a zero-mean complex Gaussian 
noise. 

Equation 0 can be represented as 

M 

yi ^ ( 6 ) 

m—1 


where Xi^km ^ ratio of complex variables, i.e. 


^l,km = 


(1 Vl,', 

hl^rn “t” ^l,m 


Let us construct 


(7) 


y = [yi,y2, ■•■, 2 / l ]^- 


( 8 ) 


C. Formal Problem Statement 

Let us define the active users set as 

S = {km-m=l,2,---M} (9) 

Our goal is to detect the set of active codes, i.e., S from the 
received noisy data y. 


D. Sparse Representation 

Typically, at a particular BR opportunity under a BS, the 
total number of BR terminals M K. This condition should 
be met by any random access based communication network. 
Otherwise, by using probability theory, it can be shown that 
two or more BR terminals will collide by selecting the same 
BR code with high probability. It is well known that if 
multiple users transmit same code then the BS cannot separate 
the corresponding users. To maintain M K, different 
approaches have been considered by the Third Generation 
Partnership Project M- For example, the collision resolution 
algorithms such as random backoff can limit the value of M 
El. In the following, we shall assume that any BR code 
will be used by at most one BR terminal at a particular BR 
opportunity Q. Let x G be a vector such that its i-th 
component will be 1 only if i G 5, and zero otherwise. Using 
the vector we can represent the data model in 0 as 


y = Cx — u + i!} 


where u = Qd:, and 


Q = 


CipAi.i 

Ci,2Ai,2 

C2,iA2,1 

C2,2A2,2 


Cl.iXl,! 


Cl,2Al,2 


Ci^kXi^k 

Cl,kXl,k 


( 10 ) 


( 11 ) 


with Qi = 0 if J ^ 5. Since M K hy assumption, the 
vector X is sparse. Thus, MUD can be seen as a problem of 
estimating the sparse vector x from the received data y. Now 
suppose that we have some prior knowledge about the number 


’ In practice, multiple BR terminals may collied by selecting same BR code, 
however, this event will occur with very low probability. Furthermore, we 
observed using simulations that this rare event does not affect the performance 
of the proposed algorithms significantly. 


of active BR users i.e., M. Then we can postulate that Xi has 
a probability density function 

Px{xi) = {I - e)So{xi) + e6i{xi) ( 12 ) 

where e = y, (5i(x) = (5o(a; — 1) and <5o(a;) is a Dirac delta 
function. 

III. Correlation based code detector 

In this section, we apply a simple correlation based algo¬ 
rithm to detect active BR codes. In detecting active codes, we 
take real part of y in ( fTOb because imaginary part contains only 
noise and interference. In the following, we shall study some 
statistical properties of A which will be helpful for developing 
the code detection algorithm. 


A. Statistical Properties of A 

In general, the components of CFR are modelled as zero 
mean complex Gaussian random variable El, i.e, - 
CAf(0, ^). Hence, the channel power can be approximated 
by its variance. The variance of CFR, i.e. cr^ ^ may not 
same for every m = 1,2, •••M. However, the BS always 
equalizes the channel power of the active users through the 
ranging procedure. In effect, the variances of CFR of all users 
must remain in a known interval due to the ranging process. 
Consequently, we can use the value of average channel power 
as an estimate of cr^ ^ for all m. Similarly, BS can estimate 
CTg ^ and (T^ ^ by using pilot-aided synchronization procedure. 
Thereby, we shall use the following assumption. 

Assumption!: = cr2,cr2. = cr^-yi g 

According to 0, Az,m is a ratio of two complex quantities. 
By applying the concept of ifTSI and using Assumption^ 
every A/^m can be modelled using a complex random variable 
A = [A]r -f i[A]i with probability density function: 

/([A]r, [A]i) = 

_, (.3) 

TT (cr^lAp +al- 2pr[\]r(Ju(rv + 2p^[\]iau(Tv) 

where al = \l-a\'^al + al + al, al = al + al, and p = (1- 
ajcr^-fiTg. Furthermore, {Xpm}i,m are independent from each 
other. In Appendix-A, we derive the mathematical expectation 
and variance of [A]r. In the following we denote them by pr 
and cr^ respectively. 

Since are real valued, the expected value of [ufr in 
(uni with respect to x and A is 

E([Mj]r) = E ( Cj^t\Xj^(\rXl, 

Vf=i 

K 

= ^[(C,.,E([A].))p,(at, = l)] 

K 

= e Pr y ^ (14) 

e=i 
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By using a similar procedure, it can be verified that the 
covariance matrix of is 


Cov([tt]r.) = e cr^ 


l|C(l,:)||i 

0 


0 

0 


where C {£,:) denotes £-th row of C. In the present scenario 
Cjj € {+1, —1}, hence 

Cov([M]r) = (15) 


TABLE I 

Correlation based multiuser detection (CMUD) 

Input: Code matrix C, a real valued decoder matrix D G 

date vector y, threshold k and an approximation of M, i.e., Mq. 
Initialization: Set T = 0, = y,j = 0. 

Loop: for j = 1, 2 ■ ■ ■ Mq 

1. Set n = 0. 

2. Construct I = {£ G {1, 2, ■ • • if} : |[D/zO-l)]r| > k }. 

3. If I 7^ 0: Set T = T U I. 

Else: Go to Step 6. 

4. Update z^^'> = C^. 

5. Update k. 

End Loop 

6. Output: Indices of active codes T = T. 


B. Correlation based multiuser detection (CMUD) algorithm 

The CMUD algorithm is summarized in Table-ll] It starts 
with an empty set T as an initial estimate of indices of active 
codes and residual vector = y. To detect the presence of 
an active BR code, the code detector coiTelates the residual 
with a decoder matrix D S In Step 2 of j-th iteration, 

the algorithm takes every column of D and computes its cross- 
product with the residual It then construct a set I such 

that 


I = {^ e {1, 2, • • • iT} : I [D/ z^^-%\ > K } 


where k is a predefined threshold. In Step 3, the set I is added 
to the active code set i.e., T = T U I. However, if it is found 
that I is empty then we assume all active codes in y has 
been detected by the algorithm and hence terminate to Step 
6. In Step 4, the residual is updated by subtracting the 
selected active codes from The residual z^^') represents 

the part of active codes that has not been detected 

yet along with noise. In Step 5, we update /c to a new value 
and repeat Steps 1-5. The algorithm needs a decoder matrix D 
and threshold k as its input. In the following section, we shall 
demonstrate some procedures for choosing those parameters. 
Furthermore, the algorithm needs an estimate of M. We 
shall describe a simple procedure in Section- HlI-EI to obtain 
a rough estimate of the parameter. When D = C, the CMUD 
will be closely related to the reduced dimension decision- 
feedback (RDDFt) algorithm 0. However, we shall show that 
an appropriate choice of D can increase the performance of 
CMUD significantly. The computational complexity of CMUD 
is very small. At every iteration, the major complexity involves 
in computing which requires LK flops. 


C. Performance Analysis 

In this section, we develop the conditions under which the 
CMUD will successfully detect all active codes. Our perfor¬ 
mance measure is based on the probability of code detection 
error by CMUD. In particular, we define the probability of 
code detection error as 


(16) 


Let us define 


K 

a = max|^D/Cj| 
i=i 

13 = max{max|D/Cj|} 
7 = max||D^||2. 


(17) 

(18) 
(19) 


Lemma 1: Let y = Ce — u + -ff where = M 

and the random vector -d has zero mean complex Gaussian 
distribution with covariance cr^I and the random vector 
is defined as in (fTol l. Let D G be a decoder matrix 

such that T)J Ci = 1; V£. Set 

T = ey-rCt + 7\/2(1 -b f) \ogK^M qgI + cr^/2 (20) 

for some given i/ > 0 and Mq > M. Assume that the decoder 
matrix D satisfies the following condition: 


T -f Mo/3 < 0.5. 


( 21 ) 


If we choose a threshold k such that 


T + Mop < K < 1 - Mop - T (22) 

then the probability of code detection error (fTfil l by CMUD 
will be upper bounded by 

Pe < (7r(l -f v) log (23) 

Proof: See Appendix-B. ■ 


D. Decoder Design 

Lemma-[T] states that for given values of M and noise 
variances, the CMUD will detect all active codes efficiently if 
the matrix D satisfies the following two conditions: 

efXra + yT + MqP < 0.5. (24) 

D/Q = l; V£. (25) 

where we denote T = 1^2(1 + v) logKyjMoa"^ + cr|/2. 
Hence, to obtain an optimum decoder, we have to design D 
that minimizes the left side of (l24l) . In Section- Ull-Dll we 
shall demonstrate a procedure of decoder design based on 


Pe = Pr{f f S\. 
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the objective. As will be seen later, the procedure requires 
to solve a high dimensional optimization problem. Note that 
the base station needs not to redesign D at every BR request 
interval, instead it will be redesigned only when the values 
of M and noise variances change significantly. Thereby, the 
proposed scheme will be compliant with the IEEE 802.16 
standards. Nevertheless, the base station always seeks low 
complex algorithm. A low complex decoder design procedure 
will be proposed in Section- IIlI-D2l which adopts a popular 
approach called minimum mean square error (MMSE) decoder 
lITSl . The MMSE decoder will not minimize the left side of 
(I 24 I 1 directly, however, exhibits moderate number of BR code 
detection performance. 

1) Dectoder-I: Let C(f) be a matrix constructed from C 
by retaining all its columns except the f-th column. To obtain 
an optimum decoder based on Lemma-[T] we need to solve the 
following optimization problem: 

{D, ,d*,/3*,%} = arg min e/r^d + + T 7 (26) 

D,q,/3,7 

subject to, DJ = 1, for ^= 1, 2, • • • AT 

K 

-d <[^ < d, for f = 1, • • • a: 

<[C(£)]^D, < /31k-i, for f = 1, • • • K 
||D^||2 <x forf=l,2,---iT. 

The optimization is convex and can be solved efficiently by 
using a Primal-Dual algorithm lfT9ll . However, the optimization 
needs to solve for L.K + 3 number of variables. In the fol¬ 
lowing, we propose a low complex decoder design procedure. 

2) Dectoder-II: The decoder design strategy is based on the 
MMSE criterion IfTSl . In particular, for every ^ £ {1,2^ - ■ ■ K}, 
we design G that minimizes E(||i:^ —[yjr-lli) along 
with the constraint Dj = 1. Here, the expectation is with 
respect to x, A and noise vector 'd. Let us define the index set 
U := {1, 2, • • • AT}. By using (fTOl i and the fact that Dj = 1, 
we see that 

Djy = X £+ ^ D/Cjd;^ — ^ D/Qj-atj -I- (27) 
jeuv jev 

Note that the value of e in (O is small. Hence, ne¬ 
glecting the values of E (xjDjCjQ/^D^Xk) ;Vj 7 ^ k and 
E (^XjJDjQjQjDgXk) ;Vj 7 ^ k, and using a similar proce¬ 
dure of (O, we obtain that 

Ei\\xe--D£[y]r\\l) (28) 

= e ||D/C,||2 + e^E(||D/[Q,7||2) 

ieuv ieu 

-2e ^ E(DjC,[Q,]jBe) + \\B£\\lal/2 

ieuv 

<e^||D/C,||2-fe^||D/C,||2E(A2) 

feu iGU 

-2e ^ \\BjC,\\lE{X) + \\T>£\\lal/2 
j&vV 

^5T>JlCDi + WD£\\lal/2 (29) 


where 5 = e(l -|- E([A]r) - 2pr) and R = I]jGU 
we have to minimize (l29l l with respect to Dg along with the 
constraint Dj = 1. Using the procedure of Lagrangian 
multiplier, it can be verified that the optimization has a closed 
form of solution: 


idIi + lal/2)-^Ci 

C/(5R-|-Ia2/2)-iC/ 


(30) 


Note that whenever the number of active users or channel 
noise variance changes significantly, we have to compute the 
inverse term ((5R-|-ItT^/2)“^ to obtain an updated D, which is 
computationally expensive. We can reduce the computational 
complexity by using the following procedure. Since L < K, 
the matrix R can be shown to be positive definite (see Il20ll ). 
Hence, the eigenvalue decomposition of R has the form: R = 
UAU^ where U^U = I and A is a diagonal matrix whose 
{k, fc)-th entry is the fc-th eigenvalue of R. Denote v = . 

Then it can be verified that 


U{SA + Ia^/2)-^v 

SKj^j+iy'^/2 


(31) 


Since (i5A + I(T^/2) is a diagonal matrix, its inverse can be 
computed easily. Eurthermore, for a given C the matrix R 
is fixed, hence we can compute the eigenvalue decomposition 
in priori. Note that the minimization of (|29] | will minimize 
the values of a and 7 in (l24l i. However, it does not consider 
minimizing jd. 


E. Estimate of total number of active users 

The CMUD algorithm in Table-Urequires an approximation 
of the total number of active users M which is unknown in 
practice. In the following, we demonstrate an approximate 
procedure to estimate M. The energy received at BS due to 
real parts of y be (using (|6ll) 


S = [y]r[y]r = 




M 


E {Cb.. — [A^ } + [Vi]r 


(32) 


Taking expectation of Q with respect to and A, we see 


M 

E - [A £,m\r^i,km } 

m—1 

M M 

— [X£,m]rCl^km}{^e,kj — [X£j]rCi^kj} 

m—1 j—1 
j^m 


E(^^) = E E 


+ (N]i-) 



(33) 


The second term of right side represents the cross product of 
different parameters which is small compared to other terms. 
By neglecting this term, it can be verified that 

E{g) = ML- 2MLnr + MLE{{[\]rf) + Lal/2 (34) 

Hence, by setting E(C/) = [y]^[y]r in (l34l i one can obtain a 
rough estimate of M. 
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IV. Sparse signal recover based approach 

A. Basis mismatched sparse representation 

The CMUD algorithm is computationally very efficient 
but cannot detect large number of active users. Furthermore, 
CMUD requires the information of variances of channel mis¬ 
match parameters i.e., (T^,(Tg,cr^. On the other hand, MUD 
being a sparse recovery problem, can be solved by using 
sparse recovery algorithms a, EH. Although the sparse 
recovery algorithms are relatively computationally demand¬ 
ing, they have been demonstrated the capability of resolving 
larger number of active users provided that the information 
of channel mismatch parameters are available. However, if 
the information of cr^, Cg, are unavailable then we cannot 
efficiently apply the standard sparse recovery algorithms a, 
EH to the data model in (doll. To overcome the channel 
mismatch limitations, we pose the user detection as a basis 
mismatched sparse signal recovery problem and develop an 
efficient algorithm to resolve the MUD problem. To understand 
the concept of basis mismatch, rewrite the data model in (fTOl i 
as: 

y — i?=[C — Q]d; 

where we need to estimate x. Since x is sparse, y — d k 
a linear combination of few selected columns of [ C — Q ]. 
Hence, we say that y — '0 has a sparse representation on the 
basis [ C — Q ]. The basis mismatch problem arises when we 
do not have accurate knowledge about a part of the basis. In 
the present scenario we do not have actual information about 

Q 

B. Solution strategy 

In this work, we adopt a total least-squares (TLS) optimiza¬ 
tion with sparsity constraint ll22l to handle the basis mismatch 
problem and obtain an estimate of x. The sparsity constraint 
TLS algorithm aims to solve the following optimization ll22l : 

{ajT.Qr} = arg min Gp(a;,Q) (35) 

aj.Q 

where, Gp{x,Q) = ||[y],. - (C - Q)x\\l + ||Q||| -f |||a;||p 

where 0 < p < 1 and the value of ^ depends on noise level. 
The sparsity constraint in dTSl l has been imposed by the £p 
norm (0 < p < 1). The optimization problem is nonconvex and 
no effective optimization solvers is available that can guarantee 
the convergence to the global minimum of dTSl) . To handle the 
problem, the work in ll22l splits the optimization into two 
parts and develop an iterative algorithm to resolve Q and x 
in alternative fashion. In particular, given the iterate at 
iteration i > 0, the iterate is obtained by solving the 
following optimization 

= arg min Gp(a;, (36) 

With a;(*) available, is found as 

= argminGp(a:«, Q). (37) 

Q 


The procedure is repeated at every iteration until a convergence 
criterion is satisfied. Note that the second optimization has a 
closed form of solution ll2^ : 

Q(*+i) = (l + ||a:«||i)'' [Ca;« - [y]r]{x^^Y. (38) 

The work in ll22l sets p = 1, which results in a convex 
optimization problem in (l36l l. However, it has been shown 
in ED^ that some nonconvex relaxations of £i norm i.e., 
(0 < p < 1) are more efficient in finding sparse solutions. 
Furthermore, they take lower number of iterations compared 
to the ti based optimization. Again, in the present case 
Xi S {1,0} which has not been considered in the optimization 
(iTSl l. In this section, we first adopt a nonconvex relaxation of 
ti norm, i.e, tp norm (0 < p < 1) for the optimization in 
(|36] |. Second, we incorporate the constraint Xi G {1,0} with 
the optimization. Finally, we propose an algorithm that applied 
an iterative reweighted least-square technique to minimize a 
sequence like dlST l. We shall demonstrate by simulations that 
the optimization yields better user detection. 

The optimization problem in (l36l l is nonconvex for 0 < 
p < 1. The work in ED suggests an iterative algorithm 
called FOCUUS to solve the optimization. At }-th iteration, 
the algorithm solves a reweighted £2 minimization problem: 

= argminS'^l^(a;, (39) 

with, 

Sl=\x, Q«) = |||[y]g - (C - Q«)a;||2 + 

where = diag{|S^'^~^^ 1^’“^, • • • , |P“^}. The 

(|39] | is a quadratic optimization problem and has closed form 
of solution. Starting from an initial estimate x^'^\ the algorithm 
update x^^'> in every iteration. The algorithm terminates when 
\x^^'> — becomes small. It has been shown in Il24ll . 

1^ that the cost function Sp'^ {x, is monotonically non¬ 
increasing on the sequence {x^^^'^i- Furthermore, the limit 
of any convergent subsequence of is a stationary 

point of the problem in (l3^ . 

To incorporate the constraint Xi G {0,1} with every 
iteration of the FOCUUS algorithm, we consider the following 
optimization problem: 

x''^^ = arg min (x, ) (40) 

Sub. to. Xt{xt — 1) = 0; t = 1,2, • • • AT. (41) 

Due to the binary constraint in (l4Tl i. the primal problem 
above is not convex. Hence, we propose a Lagrangian Dual 
Relaxation method to optimize (l40l l. The Lagrangian associate 
to (l40l)-(l4T]) is 

C{x, p) = x^V^^-^^x - (^[y]J(C - Q«) + p^)x (42) 
where, = 

+ |(C - QW)T(C - Q«) + diag{/r}^ 

(43) 

and p G is the dual variable. The dual function of (l42l) is 
g(p) = infx £(x, p). Since £{x,p) is a convex function of 
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X, we can find the minimizing x from the optimality condition 
~ which yields 

x = ^ (C(C - Q«)^[y]. + m). (44) 

Putting the value of x in (l42l i. the dual function is 


ff(M) =4 


(e(C-Q«)T[y], + M 


(45) 


Now we need to find /r that maximize (l45l l. Once we find 
fi, we can compute x using (l44l i. We apply a gradient ascent 
method to maximize (l45T l. The method requires to compute 
the gradient of ( l45l l with respect to /x. At first note that for 
any vectors r and s independent from fi, we have 

d[r^ s] 





where 




(47) 


By applying chain rule and using ( l46l l. it can be verified that 
the first order derivative of ( |45] | with respect to /x is 


dfj. 


where, 



z = {^(C - Q*)^[y]. - y); ^ = (e(C - Q.V[y]r + y) 


and Diag(A) denotes a vector constructed from the diagonal 
components of the matrix A. 

By using the results, we are now ready to develop an itera¬ 
tive algorithm to solve the basis mismatched sparse recovery 
problem. Motivated by the objectives of dTSl l and ( |39] |, the 
algorithm aims to find the sequence that will 

satisfy: 


Gp(x«,Q«) < 

where, Gp(i;^*\ = 

||[y], _ (c - qW)aW||2 + ||qW||^ + 1 [*W]T 


(49) 


(50) 


with the constraint in (EB. Note that the main difference 
of the objective in (1491 ) with ( l35T l is that instead of £i 
norm, we incorporate a reweighted quadratic term in dill). 
Similar to ll22ll . we split (l49l) into two parts. Given the 
iterate Q0 -i) the algorithm first minimizes Gp(a;,Q0-i)) 


TABLE II 

£p-TLS Algorithm 


Initialization 

1. Set Q(0) = = CT(CC^)-lj/ and i/ e [0,1]. 

For i = 1,2, ■ ■ ■ do 

2. By using the gradient of g{fi) in ED, apply a gradient ascent 
method to find /x that maximize ^(/x) and satisfies 

3. Compute by using (44). 

4. Update by using ED- 

5. If |a:0) — a:0-l)| ^ terminate. 

End for 

Output: Project onto {0,1}. 


for with constraint (HTI) . Once is obtained, the 
algorithm solves by using (l38l l. The final £p-constrained 
TLS algorithm is given in Table-HIl In Step-2, we compute 
/X by applying a standard gradient ascent algorithm m. 
In step-3, we compute using (l44l) . Note that Step-2 
ensures that By 

using the fact, it can be verified that the algorithm satisfies 
Gp(&^*\ < Gp{x^^~^\ at every iteration. 


C. Performance analysis of sparse recovery algorithms 

The optimization in dTSl l is nonconvex, hence it is difficult 
to analyse the user detection performance of the sparse TLS 
algorithm. However, if the information of variances of channel 
mismatch parameters i.e., cr^, rig, are available then we can 
represent the data model in dJO]) as 


y = Cx + w 


(51) 


where w = —u -|- £}. By using the proof of Lemma-[3] it can 
be shown that [u;]r has Gaussian distribution with covariance 
matrix crjl where cr^ 
to recover x: 


= Ma'^ + a'^. Then we can apply Lasso 


a;* = argmini||[y]p - Ca :||2 + 7 l|a;||i. 

X 2 5 


(52) 


Define the set L = {* G {1,2,---A} : |a;,i| > 0}. 
The following lemma obtained from ll2Tll describes the user 
detection performance of Lasso. 

Lemma 2: Suppose that M < L and the regularization 
parameter £ satisfies £ < , Then the Lasso 

V2<^J,Llog(A) 

will produce a unique solution a:* satisfying L C 5 with 
probability greater than 1 — 4exp(—C 2 L^/^^) where ci, C 2 are 
some constants whose values depend on the matrix C. 


V. Simulation & Results 

To evaluate the performance of the proposed algorithms, 
we develop a Monte Carlo simulation model using MATLAB. 
The simulation scenarios are on a WiMAX network based on 
the OFDMA/TDD profile. The key simulation parameters are 
outlined in Table 1. Total L = 144 subcarriers are reserved 
for the BR purpose. The modulation pulse is a root-raised- 
cosine function with a roll-off 0.22 and duration lOTg. The 
BR terminals follow a mixed channel model specified by 


















TABLE III 

Simulation Parameters 


Parameters Notation Values 


Base Frequency 

fsys 

2.3 GHz 

Channel Bandwidth 

f 

5 MHz 

Sampling Frequency 

fn 

5.6 MHz 

EFT Size 

N 

512 

Subcarrier Spacing 

A/ 

10.94 KHz 

OFDM symbol time 

Ts 

91.4 fis 

Cyclic Prefix (CP) Time 

Top 

Ts/8 

No. of CP Samples 

^cp 

64 

Frame Duration 

Tf 

5 ms 


ITU IMT-2000 standards: Ped-A, Ped-B, and Veh-A. The BR 
terminals select the channel models with equal probability. 
The mobile speed varies in the interval [0, 5] m/s for Ped-A, 
Ped-B channels, and [5,20] m/s for Veh-A. We assume that 
the variances of channel mismatch noise e and rj are equal 
i.e., (Tg = (T^. We follow the direction of ifTSll and set a = e*® 
where 0 = 5°. Four different methods are applied for code 
detection: i) CMUD with Decoder-I (Section- lIII-Dlb will be 
denoted by D-I. ii) CMUD with Decoder-II (Section- llII-D2b 
will be denoted by D-II. iii) Lasso as recommended in ll2Tll . 
and iv) (with p = 0.5) constrained TLS (Section- lIV-Bb will 
be denoted by TLS. 



Fig. 1. Performance of different algorithms with low noise variances. We 
set fje = 0.01, = 0.01 and K = 256. 

At first we evaluate the code detection performance by 
different algorithms at low noisy environment in Figure-[T1 This 
simulation will illustrate the capability of detecting largest 
number of active users by different algorithms. As can be seen, 
the D-II shows worst user detection performance whereas D-I 
performs much better than D-II. For example, with < 0.1, 
the D-II can detect only 3 users whereas D-I is capable of 
detecting upto 10 users. As being claimed. Lasso and TLS 
performs much better than correlation based approach. The 
TLS can detect 52 users with = 0.1 which is almost 5 
times larger than D-I. Finally the hgure also reveal the fact 
that, at low noisy environment, both Lasso and TLS performs 
similarly. 

We now investigate performances of different algorithms 
at high noisy environment. Figure-12] illustrates performances 


of different algorithms with measurement noise variance a^. 
Figure-l^a) compares performance of CMUD decoders. The 
performance of a particular algorithm changes slowly with 
as- For example, with M = 4 users the Pg of D-I are 0.1 
and 0.18 for as = 0.2 and 0.5 respectively. Again Pg are 
0.57 and 0.61 for as = 0.2 and 0.5 respectively with M = 6 
users. As before, the performance of D-II is worse compared 
to other algorithms. With Pg < 0.12, the D-II can detect only 
2 users. Furthermore, comparing Figure-|2a) with Figure-[T] 
we see that the user detection performances of all algorithms 
drop with increasing noise variances. Figure-jjfb) compares 
performances between Lasso and TLS. As can be seen the 
code detection performance of both algorithms decreases with 
increasing as- For TLS with M = 12, the Pg are 0.12 and 
0.24 for as = 0.2 and 0.5 respectively. The code detection 
performance gap increases with increasing the number of total 
active users. For example, with M = 20 the Pg are 0.66 and 
0.92 for as = 0.2 and 0.5 respectively. The performance of 
Lasso also changes with as- With M = 8 the Pg are 0.1 
and 0.12 for as — 0.2 and 0.5 respectively. Note that while 
Figure-[T] shows that both Lasso and TLS performs similarly 
al low noisy environment, their performance do not remain 
similar at high noisy environment in Figure-l2]]b). 

Figure-13 shows the effect of channel mismatch on the 
performance of code detection by different algorithms. Figure- 
13 a) shows that the performance of CMUD decreases rapidly 
with increasing cr^. For D-II with M = 3, the value of Pg are 
0.2 0.45 and 0.79 for cr^ = 0.1,0.15 and 0.2 respectively. 
The performance of D-I is relatively better than D-II. For 
example with M = 3 the value of Pg are 0.02 0.025 and 
0.23 for a^f = 0.1, 0.15 and 0.2 respectively. The performance 
of D-I degrades rapidly at tr^ = 0.2. In contrast, the TLS 
performs much better than CMUD decoders. In Figure-13c) 
we see that with M = 12 the value of Pg = 0.07,0.12 and 
0.54 receptively for cr^ = 0.1,0.15 and 0.2. In a similar setting 
(i.e., with M = 12), the value of Pg = 1 for D-I (see Figure- 
13a)). Figure-13]b) also exhibits that TLS performs better than 
Lasso. With M = 12 the value of Pg for Lasso are 0.09, 0.24 
and 0.74 receptively for cr,, = 0.1, 0.15 and 0.2. 

In Section- Hl-DI we have assumed that any BR code will 
be used by at most one BR subscriber and hence Xi G {0,1}. 
This assumption will be valid with high probability when the 
number of active BR users is much smaller than the number of 
available BR codes i.e., M <C iT. According to the simulation 
setting in Figure-13]b) the TLS can detect almost M — 20 users 
with high probability when ag = 0.1, cr,, = 0.1, cr,? = 0.2 
and K = 256. As the value of M becomes larger we 
should increase the value of K to remain consistent with 
the above assumption. To this aim, we investigate the user 
detection performance of D-I and TLS as a function of K in 
Figure-lH As can be seen, the TLS can detect 18 users with 
Pg = 0.11,0.12,0.2 and 0.23 for K = 200,250,300 and 400 
respectively. 

VI. Conclusion 

Multiuser detection in random access channels plays an 
important role in OFDMA wireless network to support ser¬ 
vices with bursty traffic. However, the existing algorithms 
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(a) (b) 

Fig. 3. Code detetction performance at diffemt channel mismatch noise levels erg and tr,, where erg = cr^. We set = 0.2 and K = 256. 




(a) (b) 

Fig. 4. Code detetction performance for diffemt values of K. (a) CMUD with D-I (b) TLS. We set o-g = 0.1, Un = 0.1, = 0.2. 
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suffer from low efficiency when the base station does not have 
accurate knowledge about the channel frequency response of 
active users. In this work, two different solutions have been 
proposed to alleviate the difficulty. The first solution applied 
a very simple correlation based algorithm. The theoretical 
analysis of the algorithm provides a direction to achieve 
optimum performance from the algorithm. The second solution 
applied a variant of sparse recovery algorithm which exhibits 
better user detection performance, however, with additional 
computational complexity. 


Appendix A 

Mathematical Variance oe A 
Denote the covariance matrix of [[A]r, [A]i]^ as 


V([A].,[A]i) 


^ri 
(Jri ^ i 


(53) 


where = E([A]^) - (E([A]r-))^, Ori = E([A]r[A]i) - 
E([A]r-)E([A]2). Then E([A]r) can be calculated as 

/ OO nOO 

/ [A],/([A],, [A]i)rf[A].d[A], (54) 

-OO j —OO 


bound for t, we first find a threshold T for ([A]r, [A]i) such 
that both |[A]r| and |[A]i| will remain below T with high 
probability. Once we obtain a reliable T, we can compute a 
bound on t by using (l55T l. To approximate the value of T, we 
utilize a concept called percentile level ll26l . The g percentile 
level denoted by Tg is defined as: 

/V /([A]r, [A],)d[A].d[A], > P ( 60 ) 

J-Tg J-Te 

The choice of g is dictated by two conflicting requirements: 
If g is close to 1, the estimate of threshold Tg is reliable but 
the sample space of ([A]r, [A]j) is large; if g is small, the 
sample space is reduced but the estimate is less reliable. In 
our experiments, we set g = 0.95. 

Appendix B 
Proof of Lemma-[I] 

Our proof is closely based on the works ll27l Thm. 6] and 
ll28l Thm. 4]. Let us define the event 


We evaluate the integral in polar coordinate. Let us define 

7 = 




r\l j. ■ n Pi^u r,, „ 

[AJi=fsin6'-, [AJr = fcost/+ 


(55) 


^ PrO^u n Pi^u 

a = - = - 


r = iLlkM. (56) 


(7y (7y '^^y 

Then by combining (I54ll-(|56]) along with the pdf in ( fOl ) we 
get 

,2^ roo / ^2^03 0 


E([A],) = r 
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= r 

r7rd 


'o V(^"+7')" 

dO 


dtdO 


p2-K 

(p 

cos 0 - 1 - ;7 

L 

[ V47y 

277 


T 


Next we compute E([A]^), 

/ OO ^OO 

/ [A]?/([A]., [A]0d[A].d[A]. 

-OO J —OO 

Using the similar conversion in (l55T l- (l56T l. we have 


(57) 


(58) 


mVr) 

p2t^ poo 


= r 


= r 


cos^ 6 2aT cos 6 




(f2+72)2 (t2+^2)2 (i2+^2)2 


p2t: poo 


cos^ 0 

0 


dtdO + 


LTra^ 


dtdO 

(59) 


There is no closed form solution of the first integral term 
in (|59] |. In fact, the solution become unbounded when the 
upper limiting value of f —oo. However, the solution can be 
bounded by some limiting value of f < oo. In practice, an is 
bounded and ah > ae,arj, hence there is very low probability 
of ([A]r, [A]i) taking very large magnitude. To obtain a reliable 


We shall demonstrate that S occurs with high probability and 
that under (|2TI) whenever E occurs, the active codes will be 
detected correctly. At first we provide few lemmas which will 
be used to prove the result. 

Lemma 3: Suppose that i? is a zero mean complex Gaussian 
random vector with covariance matrix a‘^1 and tt is a random 
vector as defined in (fTOl) . If (7r(l + i^) log < 1 

for some ly > 0, then the event S in dhTI) occurs with 
probability at least 1 — (7r(l + ly) logK)~^^'^K~''. 

Proof: Note that [D^itJ^is a linear combination of random 
variables {[Aj^fe],-} and by using the properties of 
discussed in Section- HlI-AI each variable [Xj,k]r is independent 
from other with mean pr and variance a^ By central limit 
theorem and using (fT4l i- (fT5T l we see that [Dju]r has a Gaus¬ 
sian distribution with mean epr variance 

MtT^||Df II 2 . Furthermore, since d has zero mean i.i.d. Gaus¬ 
sian distribution and entries of d are independent from u, we 
get that [D^ {—u + 'i?)]r has Gaussian distribution with mean 
-epr Yp=i D/Cj and variance Mal\\Tii\\l + ||D^|| 2 cr^/ 2 . 

The random variables {[Dj(—it -f d)]r}f^i are jointly 
Gaussian. Hence 

Pr{X^) = Pr{ max |[D7(—m- f t?)]r| < f) 

1 <^<A ' 

K 

> ]^Pr(|[D7(-'U-f t?)]r| < t) (62) 

1=1 


Denote = —iPrYp=i^J^j “ ALcr^||Df||| -f 

||Df |1|(T^/2. By using Gaussian tail bound, we have 


Pr{\[^Ji-u + d)]r\ < r) > 1 


ere 


ITT - pe 


2 £7^ 

e ^ 


(63) 


By using (fT7]i-(fT9]l we have 73 ^ < 

using the definition of r in (|2^ . we see that for any f G 
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{1,2,---K}: 






TT r - 


< (7r(l + logK)-^^^K-^^+’'\ 

(64) 


Now combining (l62l l. (l6?t . (l64l) we obtain that 

Pr(S) > (l - (7r(l + v) \ogK)-^/^K-^^+'''>Y . (65) 

Applying the inequality {1 — x)^ > 1 — Kx, valid for AT > 1 
and X < 1, we have 


Pr(S) > 1 - (7r(l + v) log (66) 


Lemma 4: Let y = Cg. — u + 'd where the random 

vector d has zero mean Gaussian distribution with covariance 
cr|I and mean and covariance of the random vector u is 
dehned as in (fT4l) and (ffSl l respectively. Let D G be 

a decoder matrix such that TiJ Cg = 1; Vf. Then under event 
S, we have 


max{|[Djy]j.|} < MjS + r 

inax{|[Djy]^|} > 1 - (M - l)/3 - r. 
j&s 


(67) 

( 68 ) 


Proof: At hrst we see that under the event S, 


max{|[Djy]^|} = max 


^DjQ + Dj([-M + r?]0 
ies 


< max 



< MP + T. 


(69) 


Similarly, under the event E 


max{|[Djy]^|} 

jGS ■> 


> 1 — max |D7C^| — max|D7([—M + i?]r)| 

£e5\j 

> 1 - (M-l)/3-r. (70) 


max ■ 

3^S 


1 + 


“ i 

ies\3 


D,' Cg + B {[-u - 




Now we are ready to proof the Lemma-[T] The proof follows by 
the principle of mathematical induction. Suppose that the event 
S occurs and D satishes (1211) . As dTlIi holds, we always hnd 
a threshold k satisfying (l22ll . In the hrst iteration of CMUD 
2 :(o) = y. The Lemma-U implies that by setting a threshold k 
as in (I 22 I 1 we can select at least one active code. Now suppose 
at {£ > l)-th iteration T consists of the unique indices of active 
codes only and = m with m < M. We have 

^ Cj-u + £) (71) 

i6S\T 


By using the similar procedure in (l69l l and (iTOl i it can be 
verihed that 

max {mJz^‘^-%\}<(M-m)p + T 

max {|[D 72 :^^~^ 7 r-|} > 1 — (M — m — \)(5 — r. 

Clearly the threshold k satishes (M — m)[3 + t < k < 1 — 
(M — m — l)/3 — T. Hence, CMUD will detect active code 
indices from {5 \ T}. Now consider the iterate p, when all 
active codes are detected i.e., T = 5. Then = —u + i}. 

Hence, as dSB occurs, |[D7 ^(p-i)]^| < t < k and CMUD 
will not detect any further code. 
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